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HIGHLIGHTS 


• MLM contains a non-physical formulation of the leaf-level boundary layer resistance. 

• A revised, physically consistent formulation of MLM is proposed. 

• Deposition estimates of S0 2 and 0 3 are only slightly reduced with the revised MLM. 

• HN0 3 deposition is reduced by as much as 38% during midday. 

• A reevaluation of estimated deposition of total nitrogen in the U.S. is recommended. 
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The Multilayer Model (MLM) has been used for many years to infer dry deposition fluxes from measured 
trace species concentrations and standard meteorological measurements for national networks in the 
U.S., including the U.S. Environmental Protection Agency’s Clean Air Status and Trends Network (CAST- 
Net). MLM utilizes a resistance analogy to calculate deposition velocities appropriate for whole vege- 
tative canopies, while employing a multilayer integration to account for vertically varying meteorology, 
canopy morphology and radiative transfer within the canopy. However, the MLM formulation, as it was 
originally presented and as it has been subsequently employed, contains a non-physical representation 
related to the leaf-level quasi-laminar boundary layer resistance that affects the calculation of the total 
canopy resistance. In this note, the non-physical representation of the canopy resistance as originally 
formulated in MLM is discussed and a revised, physically consistent, formulation is suggested as a 
replacement. The revised canopy resistance formulation reduces estimates of HN0 3 deposition velocities 
by as much as 38% during mid-day as compared to values generated by the original formulation. Inferred 
deposition velocities for S0 2 and 0 3 are not significantly altered by the change in formulation (<3%). 
Inferred deposition loadings of oxidized and total nitrogen from CASTNet data may be reduced by 10 
-20% and 5-10%, respectively, for the Eastern U. S. when employing the revised formulation of MLM as 
compared to the original formulation. 

Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons. 

org/licenses/by/3.0/). 


1. Introduction 

The Multilayer Model (MLM; Meyers et al„ 1998; Cooter and 
Schwede, 2000) has been used for many years to infer dry depo- 
sition fluxes from measured trace species concentrations and 
selected meteorological measurements for U.S. national networks 

(Clarke et al., 1997; Finkelstein et al., 2000; Schwede et al., 2011; 
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USEPA, 2013). MLM has been evaluated (Meyers et al., 1998; 
Finkelstein et al., 2000; Sickles and Shadwick, 2002) as well as 
modified (Wu et al., 2003a, b) but the version as described by Cooter 
and Schwede (2000) is currently employed to infer dry deposition 
fluxes from measurements made in the CASTNET network 
(Schwede et al., 2001; Baumgardner et al., 2002; Sickles and 
Shadwick, 2007; Bowker et al., 2011; USEPA, 2013). Conceptual 
components of the original model have also been used to simulate 
the bi-directional exchange of ammonia between vegetation and 
the atmosphere (Wu et al., 2009). 


http://dx.doi.org/10.1016/j. atmosenv.2014.03.056 
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Nomenclature 

h c 

height of the canopy (m) 

LAI(z) 

leaf area index (single-sided) of the canopy at 
height z (m 2 m~ 2 ) 

Ra 

aerodynamic resistance (s m _1 ) 

f a .soil 

subcanopy aerodynamic resistance for deposition to 
soil (s m 1 ) 

r b (z) 

leaf boundary layer resistance at height z (s m "’) 

r c(z) 

total resistance of canopy elements at height z (s m 
- 1 ) 

leaf cuticular resistance (s m 1 ) 

tcut 

O (z) 

total leaf surface resistance at height z (s m _1 ) 

^meso 

leaf mesophyll resistance (s m _1 ) 

G(z) 

leaf stomatal resistance at height z (s m _1 ) 

^soil 

soil uptake resistance (s m 1 ) 

S(z) 

general local sink term for trace species deposition 
to the canopy (s _1 ) 

S ong (z ) 

local sink term of the original MLM formulation for 
trace species deposition to the canopy (s _1 ) 

S rev (z) 

local sink term of the revised MLM formulation for 
trace species deposition to the canopy (s _1 ) 

v d 

whole canopy deposition velocity (m s _1 ) 

Z 

vertical coordinate from ground surface upward (m) 


MLM utilizes the ubiquitous Ohm’s Law (i.e., resistance) analogy 
to calculate deposition velocities appropriate for whole vegetative 
canopies, while employing a multilayer integration within the 
canopy to account for vertically varying meteorology, canopy 
morphology and radiative transfer. However, close inspection of the 
MLM formulation, as it is presented by Meyers et al. (1998) and as it 
is typically employed (Finkelstein et al., 2000; Finkelstein, 2001; 
Baumgardner et al., 2002; Sickles and Shadwick, 2007; Schwede 
et al., 2011), reveals a representation of the leaf-level quasi- 
laminar boundary layer resistance, r b , that is non-physical. In this 
technical note, we explore the impact of this non-physical formu- 
lation and suggest a revised formulation that is more consistent 
with the accepted conceptual model of multilayer resistance-based 
dry deposition. 


In this formulation, the aerodynamic resistance, R a , is in series 
with the whole canopy and ground surface resistances, which are 
represented by the parenthetical term in Eq. (1). At each level, the 
total canopy resistance, r c (z), is calculated as the sum of parallel 
resistances due to diffusion into the leaf stomata and absorption 
into the plant mesophyll and to adsorption onto the leaf cuticle 
(the factor of 2 appearing because adsorption may occur on either 
side of the leaf, while stomata are typically found only on the 
underside of leaves for hypostomatous plants). As represented in 
Eq. (3), the leaf boundary layer resistance r b , exists in series with 
each of these possible pathways for trace species deposition to a 
leaf. However, this representation of r b as occurring in series 
individually with each leaf deposition pathway is inherently non- 
physical. 

As a trace species molecule approaches a leaf it must first diffuse 
through any existing laminar boundary layer next to the leaf sur- 
face. At that juncture, the molecule can then either be deposited 
onto the leafs cuticle (on the top or bottom of the leaf) or it can 
diffuse into a stoma and be absorbed into the mesophyll. Given this 
conceptual picture, it is clear that the boundary layer resistance is 
in series independent of the two subsequent parallel pathways for 
trace species deposition (i.e., onto the leaf cuticle or through the 
stoma into the leaf mesophyll), so that the canopy resistance should 
be formulated as 


rc(z) = 

where, 

1 

W) ~ 

or after 

r,(z) = 

Then, 

r c (z) = 


r b(z) + n{z) 

1 2 
fs'(Z) + r me so t cut 

rearrangement, 

Lcut(f s (Z) ' r meso) 
tout + 2(r s (z) + r m eso) 

/ N tcut (r s (z) + r m eso) 
fcut + 2(r s (z) + r me so) 


and the revised local canopy sink becomes 


(4) 

(5) 

( 6 ) 
(7) 


2. Multilayer model description and revised formulation 


Within the MLM, hourly estimates of deposition velocities (V d ) 
for each deposited trace species are calculated from 


( h i 


V d = 


Ra + 


\o 


J S(z)dz 


+ - 


1 

r soi I + r a,soil 


( 1 ) 


and 

S(z) = S° rig (z) = ^ (2) 


S rev (z) 


LAI(z) 


r b( z ) ■ 


r cut (Ls(z) + r meso ) \ 
r cut + 2(r s (z) + r me so)/ 


-i 


( 8 ) 


3. Analysis of the original and revised formulations 

Consideration of limiting cases of the two formulations for r c (z) 
(Eqs. (3) and (7)) is instructive. First, in the case where the cuticular 
and mesophyllic resistances approach zero and the boundary layer 
resistance is much larger than the stomatal resistance - 

tcut —> 0 

Emeso~* 0 (9) 

r b (z) » r s (z ) 


with 


then the leaf-level resistance at height z, from Eq. (3) is given by 


r c (z) r b (z) + r s (z) + r meS o + r b (z) + r cut - ' J 1 

Parameterizations for the various resistances in Eqs. (1)— (3) are 
as specified in Meyers et al. (1998) or Cooter and Schwede (2000). 


1 1 2 1 2 

r^z) _ r b (z) + r s (z) + r meso + r b (z) + r C ut~^r b (zj + r b (z} 
3 

“ nAz)' 


( 10 ) 
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On the other hand, the leaf-level resistance from Eq. (7) 
reduces to 


rt(Z) 


( 11 ) 


Clearly, the original MLM limiting case result given by Eq. (10) is 
inconsistent with the conceptual idea that as the surface re- 
sistances become vanishingly small the overall local canopy resis- 
tance should approach the boundary layer resistance. However, the 
correct limiting case result (Eq. (11)) is produced from the revised 
formulation of Eq. (7). 

Alternatively, in the limiting case where the boundary layer 
resistance approaches zero - 

r 6 - 0 (12) 

then r c (z) from Eq. (3) becomes 

1 1 2 12 

r c (z) - r b {z) + r s (z) + r meso + r b (z) + rcut~*r s {z ) + r meso + r cut 

(13) 


and from Eq. (7) we obtain the identical limiting result 


1 

r&) 


/_n , r cu t (r s (z) + r meS o) 'N 
r C ut + 2(r s (z) + r me so )/ 

^ / tcut ± 2(r s (z) + r me so)\ _ 1 2 

\ tcut(r s (z) + r m eso) / fs(z) + r me so i"cut 

(14) 


three gaseous species, only the Vd values for HNO 3 are signifi- 
cantly affected by implementation of the revised r c formulation. 
For SO 2 and O 3 , the boundary layer resistance is smaller than the 
leaf surface resistances (r s , r cu t and r meso ) and thus the calculated 
V d values for these species are altered only slightly during the day 
(-3%). On the other hand, for HNO 3 with an effective Henry’s law 
coefficient of > 10 12 M atm -1 , the cuticular and mesophyllic re- 
sistances are essentially zero so that the deposition velocities for 
HNO 3 are reduced by as much as 38% during mid-day and 32% on 
average between 6:00 am and 6:00 pm LST. Similar results 
are obtained for the other five sites as shown in Table 1 over a 
variety of vegetation types and leaf area index (LAI) values. HNO 3 
Va reductions across the sites range from 10% at the low LAI 
site Pawnee Grasslands to 32% at Oak Ridge and 27—28% at the 
other locations. O 3 and SO 2 values are reduced by 3% or less at all 
sites. 

Interestingly, evaluations of MLM against measurements have 
not noted biases of the magnitude as are suggested here (Meyers 
et al., 1989; Cooter and Schwede, 2000; Baumgardner et al., 2002; 
Schwede et al., 2011), although Meyers et al. (1998) did report 
that mean biases for daytime estimates of Vd for HNO 3 were posi- 
tive and ranged from 0.09 cm s -1 for corn canopies to 0.47 cm s -1 
for pastures. However, the reported precisions of these biases for 
HNO 3 were much larger, ranging from ±0.88 cm s -1 for corn to 
±1.5 cm s -1 for pasture. The difficulties of obtaining accurate 
measurements of HNO 3 are well known (Huey et al., 1998; Kita 
et al., 2006), so that measurement uncertainties likely obscured 
the biases reported here. 


5. Summary, implications and recommendations 


4. Impact of the revised formulation on MLM whole canopy 
deposition velocities 

The revised formulation of r c (z) as presented in Eq. (7) was 
implemented into the version of MLM as described by Cooter and 
Schwede (2000). All other elements of MLM, including all resis- 
tance parameterizations, were left the same as specified in Cooter 
and Schwede (2000). A test meteorological data set was used to 
generate values for six monitoring sites (Table 1 ) representing a 
sampling of different vegetation species and morphology. Iden- 
tical meteorological data were used for all sites to allow a direct 
assessment of how differing site characteristics affect Vd values 
calculated by the new formulation. Fig. 1 presents a comparison of 
deposition velocities generated with the original and revised 
formulations for SO 2 , O 3 and HNO 3 for Oak Ridge, TN. For these 


Table 1 

Percent difference in 6:00 am-6:00 pm mean total deposition velocity (Va) between 
the revised and original formulations of MLM. 


Site Lat/Lon (deg) Vegetation LAI HNO 3 O 3 S0 2 

(m 2 m- 2 ) (%) (%) (%) 


Oak Ridge 

35.96N/84.28W White oak 

5.0 

-32.1 -3.3 -2.9 


Loblolly pine 



Bondville 

40.05N/88.37W Maize 

4.0 

-26.9 -2.8 -1.6 

Sequoia 

36.58N/118.75W Grass 

4.6 

-26.7 -2.4 -1.5 

National Park 

Ponderosa 




pine 



Pawnee 

40.82N/104.77W Bluegrass 

0.6 

-10.1 -1.3 -0.2 

Grasslands 




Howland 

45.17N/68.77W Spruce 

5.8 

-28.2 -3.0 -2.5 

National Forest 

White birch 



Everglades 

25.00N/80.50W Cattails 

4.0 

-28.6 -2.9 -2.0 


Sawgrass 


In this work, the non-physical representation of the canopy 
resistance as originally formulated in MLM has been noted and a 
revised, physically consistent, formulation has been suggested as a 
replacement. The revised canopy resistance formulation reduces 
estimates of HNO 3 deposition velocities by as much as 38% during 
mid-day as compared to values generated by the original formu- 
lation. Deposition velocities for SO 2 and O 3 are only slightly altered 
by the change in formulation (<3%). The MLM is used by the U.S. 
Environmental Protection Agency’s (U. S. EPA’s) Clean Air Status 
and Trends Network (CASTNet) to derive deposition velocities for 
CASTNet-based deposition estimates (Cooter and Schwede, 2000; 
U.S. EPA, 2013). Implementation of the revised, physically consis- 
tent, formulation into MLM will result in reduced estimates of dry 
deposition of gaseous HNO 3 . Dry deposition of HNO 3 represents 
35—40% of total oxidized nitrogen deposition and roughly one- 
quarter of total nitrogen (i.e., oxidized ± reduced forms) deposi- 
tion in the Eastern U.S. (Sickles and Shadwick, 2007). If the results 
of Table 1 and Fig. 1 are typical (10—30% reduction in HNO 3 depo- 
sition), estimated oxidized and total nitrogen deposition loadings 
would be reduced by as much as 10—20% and 5—10%, respectively, 
for the Eastern U. S. A re-evaluation of estimated deposition for 
total nitrogen over the U. S. should be undertaken (U.S. EPA, 2013; 
Baumgardner et al., 2002; Sickles and Shadwick, 2007) to correct 
the historical record or at minimum to provide an indication of the 
sensitivity of large-scale loadings to the way in which field mea- 
surements are interpreted. Furthermore, observations of surface- 
atmosphere exchange of reactive nitrogen species, particularly ni- 
tric acid, are needed to better constrain the parameterizations used 
to estimate total nitrogen deposition. Ideally, such efforts would 
leverage recent advances in in situ techniques (Neuman et al., 2000; 
Farmer et al., 2006; Veres et al., 2008) to directly measure nitric 
acid fluxes via eddy covariance over a range of vegetated 
environments. 
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Fig. 1 . Hourly deposition velocities for HNO 3 , SO 2 and O 3 predicted by the original and revised formulations of MLM for a test meteorological data set applied to the Oak Ridge, TN, 
forest canopy. 
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